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Statistics and thermally activated dynamics ol crack nucleation and propagation in a two-dimen- 
sional heterogeneous material containing quenched randomly distributed defects are studied theoret- 
ically. Using the generalized Griffith criterion we derive the equation of motion for the crack tip 
position accounting for dissipation, thermal noise and the random forces arising from the defects. 
We find that aggregations of defects generating long-range interaction forces (e.g., clouds of disloca- 
tions) lead to anomalously slow creep of the crack tip or even to its complete arrest. We demonstrate 
that heterogeneous materials with frozen defects contain a large number of arrested microcracks and 
that their fracture toughness is enhanced to the experimentally accessible time scales. 
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Fracture mechanisms and their relation to the mate- 
rial structure is a long-standing problem Ideal crys- 
tals are subject to fast brittle fracture - as was first 
explained by Griffith 0] - while homogeneously amor- 
phous systems exhibit slow ductile fracture controlled by 
the plastic deformation at the crack tip 0, Q . Real ma- 
terials are neither of the above. Real crystals do con- 
tain defects; but even in the ultimately disordered sub- 
stances the defects are not distributed homogeneously 
but form spatially inhomogeneous aggregates such as in- 
clusion clusters and/or dislocation pileups [3, 0,0- It is 
intuitively plausible that defect aggregates promote crack 
nucleation as the crack can settle at an energetically fa- 
vorable nucleation site. At the same time one can expect 
that random heterogeneities impede the subsequent crack 
propagation process as shown in Fig. Q] This poses the 
important question about the ultimate effect of frozen 
inhomogeneities on the fracture mechanism and, in par- 
ticular, whether frozen defects enhance or degrade the 
fracture toughness of a material. 

The inspiring work Q discussed the disorder- 
stimulated nucleation of critical cracks. In this Letter, 
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FIG. 1: Left: Sketch of an arrested crack (solid ellipse) in a 
random array of frozen dislocations (with cores represented by 
symbols _L); the dashed ellipse indicates a favorable region for 
crack nucleation. Right: Typical realization of the effective 
force —7 + Fx + fd(x), see Eq. (|2a[) . acting on the crack tip 
for short-range correlated random forces fd(x) with <5 = 3. 



we focus on the dynamics of cracks in heterogeneous 
materials and investigate fracture probabilities and the 
statistics of the fracture times. We consider both zero- 
and finite-temperature crack dynamics, the latter being 
governed by thermal activation. We restrict ourselves 
to the simplest case of cracks in a thin (quasi-) two- 
dimensional (2D) ideally elastic plate containing random 
heterogeneities. In two dimensions a crack front is a 
point - the crack tip; thus the additional effects aris- 
ing from crack front roughening are absent. Building on 
Griffith's concept of energy balance we consider crack 
tip motion governed by the dynamic energy release rate 
[U, 0] and derive the equation of motion for the crack 
tip. We include both dissipative and thermal forces, and 
the position-dependent random forces acting on the crack 
tip due to frozen material defects. We discuss three basic 
kinds of frozen inhomogeneities: (i) bond strength varia- 
tions, (ii) random impurities resulting in local compres- 
sion of the elastic medium, and (iii) frozen dislocations 
interacting with the crack. The range of elastic interac- 
tions with the crack tip increases when going from the 
type (i), to types (ii) and (iii). 

We find that the supercritical cracks can be arrested 
by heterogeneities inducing long-range elastic forces, i.e., 
by frozen dislocations (iii). We show that thermally ac- 
tivated cracks exhibit anomalously slow dynamics with 
vanishing mean velocity for all three types of disorder. 
We conclude that quenched defects effectively slow down 
the crack propagation and derive experimentally observ- 
able characteristic material properties such as the statis- 
tics of the critical stresses and the power-law distribu- 
tions of crack waiting times. This explains that in ma- 
terials containing long-range structural defects arrested 
microcracks are experimentally observable Q- 

Crack tip equation of motion. Let a single planar 
crack extend from — x/2 to +x/2 along the x direction 
of a 2D clastic medium of size L loaded in mode I by 
a uniform external stress a. In a perfectly homogeneous 
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clastic medium the Griffith crack energy is the sum of the 
clastic energy gain E e i(x) and the crack surface energy 
E s (x). The driving force for the crack tip advance is the 
release of the elastic energy quantified by the static en- 
ergy release rate. G(x) = —d x E e i = -na 2 x/2Y = Fx [![, 
where Y = E/(l — v 2 ) is the 2D Young's modulus. The 
energy release is balanced by the specific fracture energy 
7, related to the crack surface energy E s by 7 = d x E s . 
Griffith's criterion for the onset of crack growth gives 
G > G c = 7 (H, where G c is the critical energy release 
rate that can be reached by increasing the crack length 
beyond the critical value x c = j/F = 2jY '/no- 2 . 

Material heterogeneities are described by the frozen 
random forces fd{x) which are included into the Grif- 
fith's force balance. We adapt Gau ssian distributed ran- 
dom forces with zero mean value, /d(x) = 0, where the 
overbar denotes the average over disorder, and consider 
two types of forces: Short-range correlated forces (SRCF) 
fd{x)fd{x') = Asx s S a (x — x') (where a is a microscopic 
cutoff length) and long-range correlated forces (LRCF) 
fd{x)f d (x') = Asx^ 1 ^ 2 ^^- 1 ^ 2 . A s is the strength 
of the random forces proportional to the defect concen- 
tration, and the exponent 8 characterizes the elastic in- 
teraction between frozen defects and the crack. Note 
that for both types of random forces the correspond- 
ing potential energies Ed(x) of the crack tip defined by 
the relation fd(x) — —d x Ed have the same correlations 
(Ed{x + z) — E d (x)) 2 ~ Agz s+1 for z » x. It was shown 
in Ref. [f| that random bonds (random fracture tough- 
ness) result in SRCF with S = 0, impurities produce 
SRCF with 6 = 1, and dislocations induce LRCF with 
(5 = 3, i.e., 5 increases with the range of the the clastic 
interaction between the crack tip and defects ||. 

Dissipation occurs mostly near the crack tip, where 
the elastic energy transforms into heat via plastic de- 
formation [(|, and can thus be described as the local 
viscous force exerted on the tip, — r/x, with 77 being 
the tip viscosity. Including the thermal force ((t) act- 
ing on the crack tip into the force balance, we obtain 
the overdamped equation for the crack tip motion as 
rjx = G(x) — 7 + fd{x) + ((t). As we focus on slow 
crack dynamics thermal fluctuations facilitate both tran- 
sient crack healing and extension, and the system is close 
to thermal equilibrium such that it is justified to use 
(C(i)) = and correlations (C(*)C(*')) = ^r]TS(t - t') 
(ks = 1) for thermal forces. In order to complete the 
description of dynamics, the kinetic energy of the elastic 
medium has to be taken into account via the dynamic 
energy release rate G(x,x) 



G{x, x) = A(x)G(x) « (1 - x/c R ) G(x) 



(1) 



In general, A(x) decreases monotonically with increas- 
ing crack tip velocity x, starting with A(0) = 1 in the 
static limit and reaching zero A{cr) = at the Rayleigh 
wave velocity cr Q . The approximation in Eq. (JTJ com- 
plies with most experiments and will be used in what 



follows. In a homogeneous material the dynamic force 
balance G(x, x) = 7 generalizes the Griffith criterion and 
describes the energy flux into the crack tip and its sub- 
sequent conversion into crack surface energy. Including 
the viscous force — nx, the thermal force £(£), and the 
frozen random forces fd{x) into the dynamic force bal- 
ance, —rjx + G(x,x) = 7 — fd(x) — C(t), one finally arrives 
at the following equation of motion: 

V x = B(x)[-j + Fx + f d (x) + ({t)} (2a) 
B(x) = [1 + h/vcR^x/xc)}- 1 = [1 + x/b}-\ (2b) 

which is a Langevin equation with multiplicative noise 
b = tjcr/F is the characteristic crack length. For 
small cracks x <C b, B{x) k, 1, and Eq. (f2"a|) reduces to 
the usual overdamped dynamics. On the other hand, for 
large cracks, x ^ b, B(x) vanishes as B{x) ~ b/x giving 
rise to a low effective temperature. 

Crack nucleation and propagation at T — 0. In the 
absence of flaws [fd(x) = 0], cracks can be thermally 
nucleated at temperatures T > [To| , but the material is 
stable against crack formation at T = 0. However, there 
is a non-zero probability that cracks of the critical size 
x ^ Xc are "nucleated" by quenched disorder even at T — 
if the typical disorder energy gain E 2 (x c ) ~ Agxl +S , 
which scales in the same way for both SRCF and LRCF, 
compensates for the nucleation energy AE C = 7 2 /2F 
The resulting probability for disorder-induced nucleation, 
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prob[AE c + E d (x c )<0} 



(3) 



increases with the increasing disorder strength Ag. 

After its nucleation, a growing crack can get ar- 
rested by frozen disorder which, thus, can prevent frac- 
ture. Balancing the driving force and the typical ran- 
dom force developing over the distance x, (f d (x) 2 ) 1 ^ 2 ~ 
{Ag / a) l l 2 x 5 / 2 for SRCF, we find a crack arresting length 
x* = (F 2 a/ 'A^) 1 ^ 5-2 ' which characterizes the static 
force equilibrium (we focus on supercritical cracks x 3> x c 
and neglect the crack surface energy 7). For 5 < 2, 
the force equilibrium is unstable against further prop- 
agation giving rise to brittle fracture, whereas it is stable 
for 5 > 2 because sufficient driving forces become expo- 
nentially rare at x > x* , which leads to a ductile frac- 
ture mechanism. For LRCF, we find an arresting length 
x* = (F 2 / Ag) 1 ^ 6 ^ 3 ^ and the force equilibrium becomes 
stable for S > 3. 

The probability p pr0 p for the crack propagation at T = 

is given by the probability of finding a positive force on 
the crack tip for all x > x c , p P rop ~ Y\ x prob[i r x+/ ( i(x) > 

01 SRCF show two qualitatively different behaviors 
depending on the value of 5: 



Pprop 



L/a 



S>2 



(4) 
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where F a = (A^a 5 " 3 ) 1 / 2 an d F L = (A S L S ~ 2 /a) 1 / 2 are 
the characteristic forces, defined by x* = a and x* = L, 
respectively. For 6 < 2 the fracture mechanism is brittle, 
as Pprop is independent of the system size L and strongly 
increases if x* becomes comparable to the microscopic 
length a, whereas for 6 > 2, the fracture is ductile be- 
cause Pprop decreases for large L and becomes appreciable 
only if the stable crack of length x* becomes comparable 
to L. For LRCF, on the other hand, we find 



Pprop 
-,5-3^/2 



1 



-(F/F L ) 2 



(5) 



with Fl = (A^i* -3 ) which lacks an exponential de- 
pendence on L/a because the arresting forces are strongly 
correlated. The case of frozen dislocations (iii) with <5 = 3 
is marginal. Cracks are arrested for disorders A3 ;§> F 2 

f /2 

or up to high strain levels cr/Y <C fr(fd) c (fd) ~ 0-05 where 
6(fd) and C(fd) are Burgers vector and 2D concentration 
of dislocations and we used an estimate for strong dis- 
order in glass from Ref. Fracture will occur if a 
crack both nucleates and propagates with the probabil- 
ity Pfrac = .PnuciPprop- For SRCF with 5 < 2, i.e., random 
bonds (i) or random impurities (ii), the fracture proba- 
bility is finite for large L. For LRCF with 6 = 3, i.e., 
frozen dislocations (iii) , we find a complementary behav- 
ior Pnuci ~ 1 — Pprop ; which enhances fracture toughness 
because at small defect densities the nucleation is un- 
likely, whereas at large defect concentrations the propa- 
gation is suppressed. 

The fracture probability (as a function of F) equals 
the probability that F is larger than the critical force 
F c = maxrf^ — fd(x)]/x} of the given crack, i.e., 
prob[F c < F] = Pi ra c(F). For SRCF, cracks can nucle- 
ate at N ~ (L/x c ) 2 statistically independent seed loca- 
tions in the sample [H[ with critical forces drawn from 
the distribution p c (F c ) = <9FPf ra c(F c ). Then the frac- 
ture occurs at the "weakest link" if the applied load F 
exceeds the smallest critical force of all N crack nuclei. 
The average fracture force Ff rac follows from the con- 
dition 1 = Np hac (F frac ). For 5 < 2, wc find F frac ~ 
FailnN)-^ 2 - 6 ^ 2 . Because p frac (F c ) and thus p c {F c ) de- 
crease exponentially with 1/F C , see eq. (@|, the resulting 
distribution of fracture forces for S < 2 is an extreme 
value distribution of the Gumbel type, 

prob[F frac >F] ~ eM-ciNe- c ^ N ^ /2F °/ F } (6) 

with constants c\ and a. This result applies to ran- 
dom bonds (i) and random impurities (ii) a nd g eneral- 
izes previous findings for random fuse models , which 
correspond to the special case of random bonds (i) with 
5 = 0. The probability prob[Ff rac > F] in (|6]) equals the 
probability that the sample will not fracture and, thus, 
that all N statistically independent cracks are arrested. 
For LRCF, on the other hand, crack energies at differ- 
ent positions also have long-range correlations; thus ex- 
treme value statistics of fracture probabilities does not 



emerge and the fracture probability is simply given by 

Pfrac PnuclJ'prop ■ 

Dynamics of thermally activated crack propagation. 
Having established the conditions and probabilities for 
the crack arrest by heterogeneities at T = 0, wc address 
the question of to what extent these findings have to be 
modified by thermal fluctuations. While at T = 0, any 
energy barrier leads to crack arrest, thermal fluctuations 
at T > give rise to activated crack propagation. The 
equation of motion ([2a| for the crack tip resembles those 
for the overdamped motion of a particle driven over the 
one-dimensional disorder potential Ed(x) extensively re- 
viewed in Ref. [ijj]. At low temperatures, the particle 
exhibits slow dynamics due to the wide distribution of 
energy barriers giving rise to anomalously slow diffusion, 
creep, or even particle arrest fl3j | . 

We consider an ensemble of cracks arrested at T = 
by the random forces and address the question whether 
it stays arrested when the finite temperature, T > 0, is 
switched on. To this end we analyze the dynamics of a 
typical crack (since it is the typical crack that is arrested 
at T = 0, and propagating cracks represent rare events 
as follows from the previous section), which is described 
by the Fokker-Planck equation for the probability de nsity 
P(x, t) corresponding to the equation of motion (|2aj) [l4j j : 



d t P 
r/J 



-d x J 



(7) 



~Td x [B 2 (x)P] + F(x)[- 7 + Fx + f d (x)]P (8) 



After finding stationary solutions P(x) for a non-zero 
constant current J and absorbing boundary conditions 
P(L) = 0, the normalization condition L dxP(x) = 1 
determines the fracture time Tf rac = 1/J [15| . 

/>L />oo 

~ / rfa; 7TT75lrT / dze~ I{x ' z) with 
Jo TB 2 {x) J 



I{x,z) = 



rx+z 

/ du[Fu/TB(u)-A s u s /T 2 B 2 (u)] (9) 

J X 



where we took the limit of infinite L and averaged over 
disorder. The behavior of I(x, z) for large z governs the 
fracture time and is identical for both SRCF and LRCF 
141 ] . Using the asymptotics for large u, B(u) w b/u, see 
Eq. (|2b[) . we find a finite mean fracture time Tf rac for 
S < 0, whereas it diverges for 6 > 0. For 5 = 0, the mean 
fracture time is infinite for Ao > Ao, c = FTb = TrjCR, 
i.e., above the threshold disorder strength Ao jC , which is 
independent of the driving force F. 

Now wc derive the distribution of random energy bar- 
riers that govern the activated dynamics 16(. As fol- 
lows from the equation of motion (|2a)) , the effective ran- 
dom energy controlling thermal activation is <fi(x) with 
d x 4> = {—Fx — fd{x))/B(x). Therefore, we have to find 
the distribution of barriers p(E) developing in the ran- 
dom energy landscape 4>(x) of a particle located initially 
at x = Xi. This distribution evolves from the Gaussian 
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distribution of random forces /<z(x) and can be written 
as path integral in the "energy space" applying the for- 
malism that has been developed in Ref. [16[ ■ After some 
algebra, we finally find the effective barrier distribution 



p(E) 



3 -const(£/£ ) 1 -' 5/3 



(10) 



for large cracks (x ^> Xi and x ^> b) both for SRCF 
and LRCF, where E Q = b' 1 A 3 S /[3 ~ S) F~ ( - 3+5 '>/ ( - 3 -^ is the 
characteristic barrier energy. The distribution (fTQ|) at- 
tains a simple exponential form for (5 = 0. For <5 > 3 
large barriers are no longer rare, and p{E) cannot be 
normalized. Consider the crack tip starting at Xi and 
traveling over a distance x t 3> xi for the time t. For 

< S < 3, the dynamics is controlled by the high- 
est barrier, E t , it meets, obtained from the condition 

1 = {xt/a) J™dEp{E). Then E t ~ E \n 3/{3 - 5) x u and 

using the Arrhenius relation t ~ e E '/ T we find: 



x t ~ exp[(T/£;o 



,1-5/3 ^1-5/3 



't] forO<<5<3. (11) 



For (5 = 0, i.e., random bonds (i), this represents anoma- 
lously slow diffusion with the power-law dynamics x t ~ 
t T / E ° where E = A /Fb. For T < E or above the 
threshold disorder strength Ao, c , we find the vanishing 
mean velocity x t /t sa in agreement with our above re- 
sult of a diverging mean fracture time Tf rac for (5 = 0; 
(5 < leads to fast brittle fracture Xt ~ t agreeing with 
our conclusion about the corresponding finite mean frac- 
ture time. For S > the crack motion law is slower 
than any power with xt/t « 0, representing the effective 
crack arrest (meaning an infinite mean fracture time). 
For <5 > 3, the complete crack arrest, x t w 0, occurs. 
Thus an ensemble of cracks that was arrested at T = 
remains effectively arrested (in the sense of an infinite 
fracture time or zero average velocity) for heterogeneities 
with S > 0, which include frozen dislocations (iii) and 
random impurities (ii). For random bonds (i) with (5 = 
we find the anomalously slow diffusion with the power- 
law dynamics. 

Conclusion and discussion. We have derived the 
equation of motion (|2a|) for the crack tip by incorpo- 
rating effects from dissipation, thermal fluctuations, and 
frozen heterogeneities into the dynamic fracture criterion 
G(x, x) = 7. The tip equation of motion is an over- 
damped Langevin-type equation for a particle in a one- 
dimensional disordered potential. We have obtained the 
conditions and probabilities for crack arrest at zero tem- 
perature as a function of the applied stress F and the type 
of heterogeneities involved as described by the exponent 
6 [U both for the short-range (SRCF) and long-range 
(LRCF) correlated forces. For SRCF, we find complete 
crack arrest for S > 2 and the extreme value probabil- 
ity ((6]) of crack arrest for S < 2. For LRCF, cracks get 
arrested for 6 > 3, and the fracture probability is dras- 
tically reduced for hard-worked materials [171 ] containing 
frozen dislocations (iii) (S = 3). Cracks that are arrested 



at T = can propagate by thermal activation at finite 
temperatures T > 0. For heterogeneities with S > 0, 
i.e., also for random fracture toughness (i) and impuri- 
ties (ii), the disorder potential leads to the slow crack 
dynamics ([lip with zero mean velocity as the crack tip 
gets trapped in the deep potential minima. This trap- 
ping mechanism is much more efficient than the crack 
capture by crystal lattice effects [HI and explains the 
existence of arrested metastable microcracks in heteroge- 
neous materials with sizes that can be considerably larger 
than the critical crack length of the homogeneous mate- 
rial; this effect has been observed in a number of recent 
experiments 0]. Experimentally observed fracture pre- 
cursors in heterogeneous materials with power-law wait- 
ing time distributions @, EH, can also be explained 
in the framework of our theory as characteristics of the 
case <5 = of random fracture toughness (i). It remains 
an open question for future investigations whether the 
ensembles of arrested microcracks become unstable with 
respect to microcrack coalescence and slow crack growth 
by cyclic loading in fatigue experiments. 
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